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We formulate an effective theory for the atom-mediated photon-photon interactions in a two-dimensional 
"photon fluid" confined in a Fabry-Perot resonator. With the atoms modelled by a collection of anharmonic 
Lorentz oscillators, the effective interaction is evaluated to second order in the coupling constant (the anhar- 
monicity parameter). The interaction has the form of a renormalized two-dimensional delta-function potential, 
with the renormalization scale determined by the physical parameters of the system, such as density of atoms 
and the detuning of the photons relative to the resonance frequency of the atoms. For realistic values of the 
parameters, the perturbation series has to be resummed, and the effective interaction becomes independent of 
the "bare" strength of the anharmonic term. The resulting expression for the non-linear Kerr susceptibility, is 
parametrically equal to the one found earlier for a dilute gas of two-level atoms. Using our result for the effec- 
tive interaction parameter, we derive conditions for the formation of a photon fluid, both for Rydberg atoms in 
a microwave cavity and for alkali atoms in an optical cavity. 

INTRODUCTION 

Quantum physics in two-dimensions has many interesting features which give rise to effects that cannot be seen in three- 
dimensional systems. One of the most interesting two-dimensional effects is the formation of incompressible electron fluids 
that characterize the plateau states of the quantum Hall effect fH. Also high temperature superconductivity is believed to be 
essentially a two-dimensional effect. 

The interest in the physics of low-dimensional systems has motivated both theoretical and experimental searches for new kinds 
of two-dimensional many-body systems. In the case of weakly interacting Bose-condensed atomic gases, two-dimensionality 
can be reached in highly asymmetric traps [2], and quantum states similar to the quantum Hall states have been predicted for 
such systems when in rapid rotation jsj] . 

Another idea that has been advocated by one of us f^,^ is that photons also, under specific conditions in photonic traps, can 
form a two-dimensional system of weakly interacting particles with an effective mass determined by the (fixed) momentum in 
the suppressed dimension. Such a photon gas can in principle undergo phase transitions, much like a cold atomic gas, and can 
in a condensed phase sustain vortices and sound excitations, in a manner similar to that of an ordinary superfluid. 

This picture of the photons as a two-dimensional fluid has been based on the (effective) Maxwell theory of electromagnetic 
waves in a non-linear medium where only one longitudinal mode inside a cavity is excited by an incoming laser beam. The 
corresponding mean field equation has the same form as the Gross-Pitaevski equation, or the non-linear Schrodinger equation 
with a quartic non-linearity, and when coupled to an external driving field it has been referred to as the Luigiato-Lefever (LL) 
equation. The LL equation has been used to discuss the apparence of transverse patterns in the Ught trapped inside Fabry-Perot 
and ring cavities f^. 

The LL equation is a non-linear classical field equation, but it can also be interpreted as a quantum field theory with the 
electromagnetic field as an operator field. The non-linear term is then viewed as a short range ((5-function) photon-photon 
interaction. This interpretation is the basis for the photon fluid idea, and it has implications beyond the classical non-linear 
optics description. 

However, the interpretation of the non-linear field equation as a quantum theory raises several questions. One has to do with 
the dimensional reduction itself. When only the fundamental longitudinal mode is excited there is clearly an effective reduction 
of dimension, since the dynamics is restricted to the two transverse directions. This corresponds to the situation where the cavity 
is small, with a length in the longitudinal direction of the order of half a wave length. In the optical regime such a resonator is 
extremely small, and even if it can be made in principle, a simpler realization of a small resonator is in the microwave regime in 
conjunction with Rydberg atoms which can couple strongly to the microwave photons. For the cavities that are presently used in 
laser experiments the longitudinal mode is highly excited relative to the fundamental mode, and in this case two-dimensionality 
is obtained only as long as the scattering to other longitudinal modes can be neglected. 

Another question concerns the photon-photon scattering 1 7 ] . A two-dimensional delta-function interaction is only well defined 
to lowest order in perturbation theory, and in a full quantum description such a short range interaction is meaningful only as a 
renormalized interaction. This implies that the scattering amplitude is determined by a renormalization length in addition to the 
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interaction strength. In the effective photon theory this is not a free parameter, but should be determined by the full microscopic 
theory of the photons interacting with the atoms of the non-linear medium. 

In this work we will address the first question simply by assuming that only one longitudinal mode is excited. Our main 
objective will then be to examine the photon-photon interaction from a microscopic point of view. Of particular interest is 
to examine in what sense the effective interaction can be interpreted as a delta function interaction and to determine how the 
renormalized strength of the interaction depends on the physical parameters of the system. 

The approach we will take is to derive the effective photon theory from from the full quantum theory of the electromagnetic 
field and the non-linear medium, rather than from a macroscopic description of the electromagnetic field. However, we will use 
the simplified model of the atoms in the medium as a collection of Lorentz oscillators supplemented by a quartic oscillator term 
to account for the non-linearity |8]. At the quantum level, the linear Lorentz oscillator model yields "polaritons" as the coupled 
atom-photon degrees of freedom, as shown by Hopfield and others |9, id]. 

In the next section (2) we use the Feynman path-integral method to find an expression for the interaction between the physical 
modes of the coupled photon-oscillator system in terms of an effective photon action. In Section 3 we consider the effective 
theory for the low-lying transverse momentum modes in a cavity where only the lowest longitudinal mode is excited, and derive 
the corresponding two-dimensional low-energy effective action. In the following section (4), we summarize the question of 
how to correctly describe the renormalized delta function interaction in two dimensions. Then in Section 5, we relate this to 
an evaluation of quantum corrections to the effective interaction (to second order in the coupling parameter) and determine 
the leading logarithmic corrections to the scattering amplitude. In Section 6 we summarize the physical scales and discuss 
the conditions under which interesting quantum phenomena like Bose-Einstein condensation and the formation of two-photon 
bound states may take place. In section 7 we examine two possible scenarios for experimental realizations of a 2D photon 
fluid. Based on order of magnitude estimates, we discuss under what conditions a photon fluid in thermal equilibrium may form 
for millimeter-wave photons interacting with Rydberg atoms and for optical photons interacting with alkali atoms. Both cases 
might offer possibilities to observe genuine quantum effects in an interacting photon gas, although our estimates indicate severe 
constraints on the physical parameters. We close this section by discussing a third possible scenario, in which a 2D photon fluid 
forms just inside the surface of a high-Q microsphere of glass. Concluding remarks are found in Section 8. 

2 THE EFFECTIVE PHOTON ACTION 

For photons interacting off resonance with the atoms {i.e., the oscillators), the atoms have two types of effects on the scattered 
photons. There is a linear effect, where the elastic scattering off the atoms changes the dispersion of the photons. For photons 
with low transverse momentum in the cavity, this leads to a renormalization of the effective mass of the photons, which without 
the atoms is determined by the (fixed) longitudinal photon momentum hk^, as mo = Tik^ /c. The other effect is the non-linear 
or anharmonic effect of the photon-atom interaction, which gives rise to the effective photon-photon interaction. 

There are several ways to derive the effective photon action from the full quantum theory of photons and oscillators, all based 
on the assumption that the non-linear term is small and can be treated as a perturbation. One way is to solve, as the first step, 
the linear part of the problem exactly by diagonalizing the Hamiltonian. As the next step the non-linear terms can be expressed 
in terms of the transformed, decoupled variables; in this way the resonance problems, which would appear in a more direct 
perturbative treatment, are avoided. However, another simpler approach which we will adopt here, is to derive the effective 
action of the electromagnetic field by use of the Feynman path integral method, where a (non-local) field transformation yields 
the result without any matrix diagonalization. To check the result of this method we have also performed, in Appendix A, a 
decoupUng of the linear variables by matrix diagonalization, and show that this can be done in a way that is substantially simpler 
than the standard method. 

In this section we do not impose the cavity boundary conditions, which later will be used as constraints on the effective photon 
modes in order to derive a dimensionally reduced theory. 
We start from the total classical action 

S[A, R,] = [A] + Sm [R^] + Si [A, R,] (1) 

which is the sum of a free photon part Sj, a matter part Sm, for the atoms, and an interaction Sj which describes the their 
coupling to the radiation. Ri describes the atomic degrees of freedom, which here are given as the displacement vectors of a 
discrete set of oscillators labeled by i. In the following we will put h — c — 1. 
The free photon part, S^, is given by the Maxwell term 



5^ = 1 /" d^rdt[E^ ~ B^], (2) 
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where E = —A and B = W x A with A satisfying the Coulomb gauge condition V • A = 0. The oscillator part includes an 
anharmonic term, 

Sm ^ J [y^' ~ T^°^^ " ' 

and the atoms are coupled to the electromagnetic field via a dipole interaction, 

Si^ J dtY,qA{n)-R,. (4) 

i 

Here fi is the spatial position of the i:th oscillator, and we shall furthermore assume that the atoms are uniformly distributed in 
space with a number density p, so that discrete sums can be replaced by integrals over a continuous position vector r. Since 
only the transverse part of R couples to the photon field, it is consistent to neglect the longitudinal part and impose a Coulomb 
"gauge" condition also on the oscillator, i.e., V • i? = 0. 

The effective photon action Sej / [A] is defined by the following (path) integral over the oscillator variables Ri 

i 

To perform the Ri integrations, we go to Fourier space and expand to lowest order in A, 



1 - 
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■ / n ^ 2^<5(^ uj^)R,{uj,) ■ R,{cj2)M^3) ■ i?.(^4) + 0(A')] . 
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The first term is evaluated directly by completing the square and making the shift 



Ri{uj) — > Ri{ijj) — iujj7Ai'DQ{u!), with 'Do{llj) as the retarded propagator. 



Vo{uj) = — • (7) 

— + IUJ€ 



(8) 



It yields the quadratic part of the effective action, 

where in this expression we have taken the continuum limit ^ - ^ p J <Px and also introduced the effective plasma frequency 

7? = [pq^/Mf/^. 

The quartic R term in (|6} is most easily evaluated by again performing the shift in Ri{uj), to give A dependent terms of the 
form A"*, {A ■ A){R ■ R) and {A ■ R){A ■ R). The A"*" term can be directly re-exponentiated and gives a quartic contribution to 
the effective action, which in the continuum limit is 

^iff = [d'^ [U \^^n VoiuJn)] S{J2 ^n)Ai^l,x) ■ A{uj2, x) l(c.3 , ^) ' A{^i. x) . (9) 

P J n=l L -I n=l 

The terms proportional to {A ■ A){R ■ R) and {A ■ R){A ■ R) can be evaluated using the formula, 

I ]JVR.,e'I^^^^^^^--^(-'--«^^^^-^R]{^^^^^^ , (10) 

i 

where A/" is a field-independent normalization factor They give, in principle, a correction term to the quadratic action, but due 
to the integration over ^^{uji) this correction term vanishes. 

Note that the expansion and re-exponentiation of the non-linear term will generate correction terms, but these are higher order 
in A and will be neglected. Thus, to first order in A, the effective action is given by the quadratic part (|8j and the quartic term 

i). 
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The effective action defined by (|8} and (|9j corresponds to a Lagrangian that is non-local in time. However by a further 
transformation it can be brought into a local form. We first note that the quadratic part of the action defines a modified dispersion 
equation 



2 7 2 2 

uj — k — rj 



= 



(11) 



with solutions 



(12) 



This equation defines the dispersion of the "polaritons", i.e., the two decoupled degrees of freedom of the linear problem which 
mixes the photon and dipole variables. For fc^ < cJq + rf' uj^ represents essentially the photon mode and the dipole mode, 
whereas for fc^ > ciJq + -q^ the interpretation of the two modes is reversed. In the intermediate interval with /c^ w Wq the photon 
and the dipole modes are strongly mixed. 

The following field transformation is now applied 



^2 _ 



\ Cj2-t^±(fc2) 



and this gives for the quadratic part of the action 

,(2) _ 1 ( d^k 



S. 



eff 



(27r)3 



^ {k, oj)* {iu^ - u;l) {k, u;) . 



(13) 



(14) 



The dependence on cj^ shows that the transformed action corresponds to a Lagrangian that is local in time. The non-locality in 
time has been traded for a non-locality in space, but this is less problematic in a Lagrangian formulation. Note, however, the 
ambiguity in the transformations (I13> . depending of which one of the solutions lj± we choose. Clearly, the relevant choice is the 
one which fits the energy of the photons in the effective theory. This means ui- for the case of red detuning (energy below ujq) 
and UJ+ for blue detuning (energy above loq). 

When the transformed field is introduced in the quartic part of the effective action we make a further simplification by assuming 
that the fields satisfy the dispersion equation of the linear problem. This allows the following substitution 



iu!'Do{uj)A{k, uj) 



lim 



1 



^ uj'^ — k'^ — rj'^ 



E^^\k,u;) 



where E'^^'' (k, uj) — iuj±A^^'> (fc, w). For the interaction part of the action this gives the following expression, 



(15) 



S. 



(4) 
eff 



(27r)4?74A 



n 



(Pkr, dWri 



(27^)4 V(a;±(fc„)2_^2)2+^2^2 



5{Y,kn)5{Y. 



x3^\ki,uji) ■ 3^\k2,u;2)E^^Hh,uj3) ■ E^^\kA,uji) . 



(16) 



The application of the dispersion equation to the field variables of the interaction term can be justified when this term is used 
perturbatively, with the fields satisfying the field equation of the unperturbed system. However, it is interesting to note that the 
expression ( I16> in fact is valid beyond this approximation, as is demonstrated by the diagonalization of the quadratic problem 
performed in Appendix A. 



3 DIMENSIONAL REDUCTION AND THE EFFECTIVE 2D THEORY 



Due to the boundary conditions imposed by the mirrors, the component of the photon momentum normal to the mirrors (the 
longitudinal momentum) is quantized at discrete values. We assume an idealized situation with infinite flat mirrors, thus the 
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FIG. 1: Dispersion curves for the photon like and dipole like excitations. Fig la shows the curves for red detuning (A — ljo ~ m > 0), 
where the photon branch corresponds to the curve ui- for low momentum k and to ui+ for high momentum. Fig lb shows the curves for blue 
detuning (A = tJo — m < 0), where the photon branch corresponds to u!+ for all momenta. The curves are shown in dimensionless units with 
|A| — 1. The parameter values are ujo = 100 and r; = 0.1. 



longitudinal momenta are quantized as — irn/ L, with L as the distance between the mirrors and n as an integer We also 
assume photons to be fed to the cavity by a laser (or maser) tuned close to resonance with one of the modes, either slightly 
below (red detuning), or slightly above the resonance (blue detuning). However, we do not take the effect of photons entering or 
departing the cavity explicitly into account, and in this sense we consider an idealized situation with perfectly reflecting mirrors. 
All the photons inside the cavity are assumed to be trapped in the same longitudinal mode, and throughout the paper we will 
assume this to be the lowest mode n = 1. 

The transverse components of the photon momentum k we assume to be restricted to small values, kr << kL- The disper- 
sion of free (non-interacting) photons inside the Fabry-Perot resonator becomes essentially that of 2D massive, non-relativistic 
particles 01, 

I p 

oJk = \Jkl + klKkL + (17) 

with the longitudinal momentum k^ playing the role of the photon mass. The dimensional reduction is then based on the 
assumption that only one longitudinal mode (the lowest) is excited, and that scattering to other modes can be neglected. We 
should stress that this does not mean that higher modes are not important as virtual states in the perturbative expansion - in fact 
they are. 

For simplicity we shall in the following refer to the transverse momentum simply as kx = k and the longitudinal momentum 

as fci = tt/L = mg. 

When the coupling between the photons and the oscillators is taken into account, the dispersion equation is given by ( II 2t . and 
the relation between k and uj{k) is no longer so simple. In Fig. la w_(fc) and u!+{k) are shown as functions of the transverse 
momentum for red detuning, which is the case we will first consider It displays how for small momenta, (fc) corresponds to 
the "photon branch" with quadratic dependence on k, while for large momenta the photon branch is represented by lu^ (k) . In 
the following we will simply refer to excitations of this mode as "photons" and the other one as "dipoles". Due to the mixing 
there is an avoided level crossing at intermediate momenta, where there is no clear distinction between the photon and the dipole 
mode. 

For the case of interacting photons, just as for non-interacting photons, a low-momentum description can be made where 
the photons appear as non-relativistic, massive particles. Thus, when the photon frequency is separated from the resonance 
frequency of the oscillator mode by a detuning gap A = luq — niQ, and the transverse momentum is restricted by 

fc^ < 7noA < ujI , (18) 

then the previous expressions for S^^^j, and 5^^^ ((|8j and (|9jl) define a low momentum effective action for the photons, with a 
dispersion relation of a non-relativistic form 



fc2 



(19) 
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In the following we shall in addition assume weak coupling between the photons and dipoles, in the sense 

77' < . 



The effective photon mass is then given by 



77 



2 mo 
2A 



(20) 



(21) 



with a small renormalization of the mass due to the interaction with the dipole field. Note that the weak coupling condition iZOi 
is not essential for the non-relativistic description of the 2D photons, but is introduced to simplify the calculations. In physical 
realizations of the 2D photon gas one may also have to consider the case of strong mixing of the photon and dipole degrees of 
freedom, as discussed in the section on Rydberg atoms below. 

With the longitudinal momentum fixed to mo and Lu{k) approximated by m, the field variable A'-^^ can be written as. 



a=±l 



sin(moz) \e--i--i->^-r) ^^(k) ^^^k,t) + , 



(22) 



where ea{k) are the polarization vectors and (paik, t) are the corresponding field components, which now only depend on the 
transverse momentum k. Note that both the frequency and the momentum of the longitudinal mode have been extracted from 
(t>a{k, t) in order to express this as a slowly varying field. 

When the assumptions about small transverse momentum and weak coupling are imposed, the quadratic part of the effective 
action gets the form (to order fc^). 



a=±l L 



1 

2m 



(23) 



Here we have neglected terms proportional to (slowly-varying field approximation), and V is now the two-dimensional 
gradient. Note that the two polarization directions appear as two species of particles. The action has the standard form of a 
non-relativistic, free field theory. 

We will now consider the interaction term. In the same approximation as used above we have 



1 



^{uj^{ky-ujlY+r^-^uol 2m|A| 



(24) 



If this fc-independent expression is used in the interaction part of the action, (|5Jl, and the i?^ ^ field is expressed in terms of (j)a, 
it simplifies to 



/ d\dt 

2 m 



a=±l 



where g is the (bare) interaction strength given by 



9 = 



1 3 A 



(27r)2 16 p 



(!)■ 



(25) 



(26) 



The first term of \25\ can be interpreted as a delta function interaction between photons with the same helicity, the second be- 
tween photons of opposite helicities. In the simplest case, with only one type of photon polarization {(f)a — (p), the corresponding 
interaction Lagrangian simplifies to 



(4) 



2m 



(27) 



Thus, with the approximations used, we reach a form of the effective photon Lagrangian which agrees with the nonlinear 
Schrodinger equation previously derived from the classical field theory of a dimensionally reduced Maxwell field interacting 
with a non-linear medium 1 11 ]. The photons behave like 2D massive particles with a repulsive, pointlike interaction. 

In the case of blue detuning, i.e., A < 0, the situation is quite different. Typical dispersion curves are shown in Fig. lb, and 
we see that now the uj+{k) branch is photon-like both for small and large momenta, and there is no avoided level crossing, but 
only a level repulsion at small momenta. 
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However, also for blue detuning a low momentum effective action can be given where all the above formula, derived for red 
detuning, still hold for scattering between the photons. But an important change is that the photons now are the "high energy" 
particles relative to the dipoles, and it is energetically possible for them to "decay" into the low lying dipole modes. At tree 
level (to first order in the interaction) the dangerous process is 77 7D which for blue detuning can conserve both energy and 
momentum. The corresponding interaction piece in the Lagrangian can again be read from il6\ . but now with one of the fields 
as E^~'> (dipole field) and the three others as (photon field). In the low momentum approximation it is 

4t = -^^(r'^)(r0), (28) 

where cj) is the photon field and <j) is the dipole field. 

An important quantity for assessing whether a gas of blue detuned photons can be maintained in the cavity, is the ratio 
between the cross section of the above ('inelastic') decay process, and the normal ('elastic') scattering 77 77 induced by the 
interaction (I26> . A straightforward calculation yields, 

^tot(fc) - li'^afl (29) 

\ '/ / 

where k is the momentum in the center of mass. Thus, demanding crl"f- <^ crfot^ implies the condition 




^«(I)VT- 

We shall return to numerical estimates of the physical parameters in section 7. 

We conclude this section with some comments on Galilean invariance. The effective action, given by (|8} and (|9j, is neither 
Lorentz nor Galilean invariant, since the relativistic photons are coupled to dipoles defined in a fixed frame. Nevertheless the 
effective theory defined by ( I23> and (125 > is Galilean invariant due to the low momentum approximations made in the vertices. 
In the next section we will consider loop effects where the low momentum approximations are not any longer valid, and it is far 
from obvious that the resulting corrections to the effective theory will respect the Galilean invariance. As we shall see, however, 
the leading corrections do have this symmetry, so the interpretation of the photons as a nonrelativistic Bose system has validity 
beyond the Born approximation. 



4 RENORMALIZATION OF THE 5-FUNCTION INTERACTION. 

In a many-particle interpretation the interaction Lagrangian Ml\ corresponds to a delta-function potential 

^w-^E-^^'H^;,) (31) 

where Vij is the two-particle relative position. However, it is well known that a pure delta function interaction in dimensions 
higher than one is not well defined beyond first order in perturbation theory. In two dimensions the second order term gives rise to 
a logarithmic divergence in the scattering amplitude. To make the delta function interaction meaningful requires regularization 
of the interaction and renormalization of the interaction strength. As discussed in ref.|12] the form of the s-wave scattering 
amplitude for such a renormalized interaction is 



f(k) = -(27r)^/-^ -^1- ^ (32) 



1 



27rfc 



0^ 



where g,. is the renormalized coupling constant and fir is a new parameter that is introduced by the renormalization (the renor- 
malization scale). The corresponding phase shift for small is given by 5o = ngr/i, which is the Born approximation value 
when g = gr- For small valus of k it approaches the universal expression, (5o ~ vr/ In/c^, that is common for a large class of 
short range potentials 

Formally, the delta function interaction in two dimensions is dimensionless, i.e., it scales as the kinetic energy. However, 
the renormalization breaks the scaling symmetry and introduces a length scale through the parameter /i^. This is similar to the 
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situation in QCD, where the effect is referred to as dimensional transmutation lfl4ll . One should note that the two parameters 
Qr and fir are not independent. Thus, may be fixed as the bare parameter g and all the effect of the renormalization may be 
absorbed in n^, or may be viewed as depending on /i^, where /i^ is chosen to match the physical momentum interval. In the 
latter case is referred to as an effective (or running) coupling constant. The explicit dependence on /i,. is given by 

- =i-ln4 (33) 
5r 9 Mo 

with /^o as a constant. From this expression we notice that for large momenta (large ^r) the effective coupling constant goes to 
zero, so this is a quantum mechanical analogue of the asymptotic freedom of QCD. Note also the curious fact that for sufficiently 
large k the effective interaction is always attractive, irrespectively of the sign of the bare parameter g, whereas it in the other 
limit is repulsive. 

The above discussion refers to a situation where the (f)'^ theory is treated as a fundamental theory where gr and fir are free 
parameters, to be determined by experiment. However, treated as a an effective (low energy) theory they are in principle 
determined by the physical parameters of the complete system. For example, in the case of a quasi two-dimensional atomic Bose 
gas in a highly asymmetric trap, the renormalization scale of the two-dimensional theory is essentially given by the extension of 
the trap perpendicular to the plane in which the atoms move I^Tsi flill . 

In the present case the renormalized interaction strength may be determined by taking more explicitly into account the ef- 
fect of the dipole degrees of freedom. This we do by using Schrodinger perturbation thery, with the interaction Hamiltonian 
extracted from the the action il6\ . to calculate contributions to the scattering amplitude beyond the Born approximation. With 
the expression for the scattering amplitude given by i32\ . which we assume to be correct for low momenta k, we note that the 
renormahzation scale can be determined from the contribution to second order in g (with — g)- In such a second order 
calculation of the scattering amplitude the contributions from the dipole mode cannot be neglected, since the intermediate states 
are not restricted to low momenta. Thus, both the field modes E'-^'> are included in the calculation, and the exact expressions for 
the mode frequences uj±{k) are used rather than the low momentum approximations. 

In the following we present a simple calculation of the leading contributions to the scattering amplitude. The result shows that 
the form of the amplitude is as expected and it gives an estimate of the renormalization scale /i^. In Appendix B we perform a 
more complete calculation of some of the non-leading terms in the scattering amplitudes and give the corresponding expressions 
for the scattering amplitude, for both red and blue detuning. 



5 THE SCATTERING AMPLITUDE TO 0{g'^). 

Consider the T-matrix, related to the scattering matrix by Sfi = Sfi — 2TriS{Ef — Ei)Tfi, where Ef and Ei are the energies 
of the final and initial states. To second order Tfi has the form 

Tf, = {k^k^\H„,t\kik2) + - / rf Pirf P2 . , jt—. 1 \ , ■ (34) 

where we here have simplified the notation by the sum over field modes, longitudinal momenta and polarization variables in the 
intermediate state. The form of the interaction matrix element is 

1 Ary'' 1 

{PlP2\ Hint \klk2) = (27r)3 2pA'^ ^^^^ + k2- Pi - P2)((5riiri2 + 2^nilSn2l - <^nin2-2) 
V i^iPl) - ^o) + V^O i^{P2r - ^o) + ^ '^0 

where we have approximated the energy of the incoming photons (fci and k2) by m, since these are in the lowest longitudinal 
mode. These photons have the same polarization vector e, while the particles ("polaritons") in the intermediate state have 
polarization vectors ea{pi) and ei}{p2)- The quantum numbers ni and n2 determine the longitudinal momenta of the particles 
in the intermediate state. 

In the low-momentum approximation the T-matrix element is independent of the sum of the transverse momenta, K = ki+k2, 
i.e., it is Galilean invariant. This follows since a;(fci) and uj{k2) can be approximated by m in all places, except in the energy 
denominator when the intermediate particles are also in the lowest longitudinal mode. In that case the pole at P1+P2 = kf + fc| 
makes the fc-dependence important. However, due to momentum conservation, we have in the low-energy approximation 

uj{ki) + Lj{k2) - uj{pi) ~ uj{p2) = -^{kf + kl - pj - pI) 
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FIG. 2: Second order contributions to tiie scattering amplitude. Only contributions from diagram (a) are included in this paper, since the 
contributions from (b) and (c) are suppressed by a factor A/m. 



where k = ^{ki — fc2) is the relative momentum. The expression shows that K can be absorbed in the integration variable pi. 
Thus, the T-matrix has the following momentum dependence 

Tf^ = 5{Kf-K^)T{kf,h) (37) 

where Kf and Ki are the sum of momenta for the outgoing and incoming particles and kf and ki are the relative momenta. For 
pure s-wave scattering (delta function interaction) the reduced T-matrix element T only depends on the magnitude of the relative 
momentum, 

T{kj,ki)=T{k') (38) 
and is related to the s-wave scattering amplitude through 



/(fc^) = -(27r)^my^r(fc2) (39) 

Since we are interested in the behaviour of the scattering amplitude for small transverse momenta, in the first order expression 
we simply put them equal to zero. With all photons in the same helicity state, the first order contribution is 

^i(fc') = TT- (40) 

in accordance with the expression for the low-energy Lagrangian ( I27> . To second order there are three diagrams shown in fig. 2. 
Diagram 2a corresponds to two particles in the intermediate states, while diagrams 2b and 2c correspond to four and six particles 
in the intermediate states. Thus, the contributions from diagrams 2b and 2c are suppressed by the energy denominator, since the 
energy difference between the two initial and the four or six intermediate particles necessarily has to be large on the scale set by 
the transverse momentum. For this reason we shall only consider contributions from diagram 2a. Note that there are two types 
of particles in the intermediate state, characterized by energies lo^ and lo^. As an important point also note that the transverse 
momenta of the intermediate particle states cannot be assumed to be small, and also excitations to higher longitudinal momenta 
have to be considered. We now discuss how to calculate the leading contributions to diagram 2a. 

At low momenta, there is a potentially large contribution to diagram 2a, when the energy denominator vanishes, and as 
expected this will give rise to the logarithmically infrared divergent term ~ In fc^ in ( I32> . This term is dominant in the limit of 
asymptotically small momenta where the approximations leading to J16t become exact. 

The importance of the high momentum contribution comes from the fact that w_ does not increase with momentum, but rather 
approaches the resonance value ojq, c.f. Fig. 1. Thus, even if intermediate states with frequencies close to cjq are considered 
as highly excited relative to the low energy photons, the large number of dipoles may make contributions from these modes 
important. In fact, if the dipoles are treated as a continuum, the integral over intermediate momenta will diverge. In reality we 
know that there is a physical cutoff related to the discreteness of the system of dipoles. We introduce this simply as a cutoff in 
momentum at a value corresponding to the (average) distance between the dipoles. 

With a clear separation of the scales in the momentum integrals the leading contributions from high and low momenta can be 
estimated separately. In order to see how this works we examine the following toy problem. Consider the integral 

1 + x"^ 

dx—-^ (41) 
which can be evaluated exactly, to give 

I = ^' + log [xlrr^^] + P' log [^^] - log [-r^-u] ) (42) 
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However, assuming 

/ < < < (43) 
the integral can be estimated by the following approximation. 







dx —+ / dx — (44) 

X + a; + j„j2 a; + a; + A^ 



~ / "a; ^ + / da; 1 

4 2 

- -A^l^A^ + ^ ' 

where the low and high momentum contributions are treated separately. This expression reproduces the exact result up to 
©(m^/A^). Below we will examine the leading contributions to the second order scattering amplitude in this way, by evaluating 
separately the contributions from low and high momenta. In this calculation the detuning parameter and the photon mass will 
play the role of A and m in the toy problem. We refer to Appendix B for a more complete treatment. 



The high-momentum contribution 

For large momenta the important contribution comes from the term with two dipole excitations in the intermediate state. For 
these excitations we have have uj{p)'^ ~ oJq and the momentum integral is divergent without the cutoff. The only effect of the 
coupling between the photons and the oscillators appears in the denominators of the form 



1 1 



(45) 



where it prevents the expression from diverging when io{p) Wg- 

We can neglect contributions from the transverse momenta of the scattered (external) photons, since these are smaller by 
0[{k/mY] relative to the leading term. This means that the contribution gives rise to a fc-independent renormalization of the 
interaction strength. We also neglect terms that are higher order in the coupling strength -q. With Loik) — m for the external 
photon states and a;(p) = o^o for the intermediate states the energy denominator is approximated by 

(46) 



u){ki) + Lj{k2) - uj{pi) - uj{p2) 4mA 
and the matrix elements of the interaction get a simple form 

Krf' 

{pi,P2\H,r,t\kik2) = — e*-ea(pi) e * • e>(p2) 

X S{pi+P2-ki-k2)(Sn,n2+^SnilSn2l-^ni{n2-2)) (47) 

For high momenta the summation over the polarization vectors gives simply a factor 3/4 for each intermediate particle (as 
discussed in Appendix B) and the momentum integral and sum therefore gets trivial, with a momentum-independent matrix 
element. Integrating over transverse momenta and summing over longitudinal momenta gives 

T2a{k') = -^^T,^aT,, (48) 
64 m^A 

where a is a new dimensionless parameter. Since the dimensional reduction is not effective at high momenta (we have to sum 
also over the longitudinal momenta) this parameter is a characteristic of the full three-dimensional theory, and is a measure of 
the importance of renormalization of the nonlinear effects for intermediate momenta, i.e., for g In fc^ A^ <C a. 

In order to obtain the finite result ( I48> we have introduced a cutoff, Pmax — t^/Iosc, in the momentum integration, where 
lose = P^^ is the distance between the oscillators, and a cutoff in the discrete longitudinal momentum variable at — L/losc ~ 
Lnpmax- This means that we have [pmax/'^)'^ = p, with p as the 3D oscillator density. 
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The low-momentum contribution 

The important low momentum contribution comes from the term with two photons in the intermediate state. The energy 
denominator then vanishes when the momenta of the intermediate photons match the ones of the external photons. Since 
all photons then are low momentum photons, for the leading contribution we can replace oj{p) by the low energy expression 
m + /2m (and ijj{k) by m + k"^ /2m) to get the energy denominator on the form 

1 2m 



w(fci) + Ld{k2) - io{pi) - uj{p2) kf + k^-pj- pI 



(49) 



Here pi and p2 are the transverse photon momenta of the intermediate state. In this case only the lowest longitudinal mode has 
to be included in the intermediate state. In the same way as for high momenta there are corrections, but they are supressed by 
factors if' /w? or A^/m^, and we shall neglect them. 

We note that an integration of the intermediate state momentum of this term alone gives rise to a logarithmic ultraviolet 
divergence. Eventually this divergence is of course cut off by the interparticle distance, but before that the integrand is suppressed 
by the factor 



{lo{pY - ujIY + Tf^l y {2mA - p^y + r]^u;^ 



(50) 



which will introduce an effective cutoff in the momentum integral at 2m| A| and thus provide a scale for the logarithm. If this is 
introduced as an explicit cutoff, the momentum integral, after the angular integration, gets the form 



2m\A\ 



1 1 



I = 2tT I dp^—^ t: = — — r ( In — — - — ITT | (51) 

Jo ^{k^~p^ + ie) (2mA)4 V 2m|A| / ^ ^ 

There will also be a constant (k independent) contribution, but as shown explicitly in Appendix B this is generally small com- 
pared to the leading high momentum contribution ( I48t . With the relevant constants and symmetry factors included, the logarith- 
mic contribution to T is, for red and blue detuning. 

There is also a term corresponding to one photon and one dipole in the intermediate state, but the real part of this is subleading 
relative to the terms already included and can therefore be omitted. However, for blue detuning it has an imaginary part 

r2^J-(fc2) = j^g(|.)Vi. (53) 

Although small compared to the leading contribution, this is the dominant imaginary part corresponding to the decay process 
77 7Z? discussed earlier. As a check on our calculations, we have verified that this imaginary part of the scattering amplitude 
is correctly related to the total inelastic cross section in ( I29> . 

The scattering amplitude 

Combining (I40> . (I48> and ( I52t we get the following approximation for the (one loop) scattering amplitude corresponding to 
red detuning 



2mn 



2mA 



(54) 



The expression is consistent with the expression for the scattering amplitude of a renormaUzed delta function interaction, when 
expanded to second order in the coupling strength. Resumming the diagram 2a as a geometrical series in fact gives the full 
scattering amplitude (I32> . if we set gr — g and define the renormalization scale /i^ = no by 



Mo 



I = 2m|A|e-T^ . (55) 
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The same expression is valid for blue detuning if we neglect the effect of scattering 77 jD. Note, however, that the sign of 
a is different in the two cases. The value of the exponent is 

-l£ ,56) 

and we note that this (in absolute value) will be much larger than 1 with the assumptions about parameter values which we have 
made. 

The choice of renormalization scale, /io in ( I55> is however misleading in that the logarithm becomes large. The definition 
/ir = 2m| A| amounts to a more natural choice since the interpretation of the photons as massive non-relativistic particles only 
makes sense for k"^ /{2m) < | A|. The corresponding value for the renormalized coupling constant is 

9r = (57) 
1 — a 



and relative to the bare (first order) coupling constant g, the change in gr remains small as long as the logarithm is small and 
a A/(r7i^A) << 1. But depending on 1 
renormalization effect. For large a we have 



a A/(r7i^A) << 1. But depending on the parameter values, a may in reality become large and give rise to a significant 



g 1 m^ri'' 

9r^ ^^^^ (58) 

and we note that in this limit the effective interaction is independent both of momentum and of the anharmonicity parameter A 
of the oscillator spectrum. The detuning parameter A (and not the anharmonicity parameter A) now determines the sign of the 
effective coupling, with repulsive interaction for red detuning and attractive interaction for blue detuning. 



Comparison with the Kerr nonUnear susceptibiUty coefficient x'^^ for two-level atoms 

Since the renormalization parameter a is proportional to 1/A, it easily gets large for small detuning, A << wo> as shown 
explicitly for the case of photons interacting with Rydberg atoms, in the section below. The renormalized coupling constant gr 
(which then is much smaller than the bare coupling constant g) should then be interpreted as the physical interaction parameter. 
One should, however, note that the expression we have found for gr is not based on a systematic expansion in A, but rather 
by resumming "dangerous" terms in the expansion. There will be other contributions to the renormalized coupling, but these 
are parametrically small, i.e., supressed by powers of small ratios like rj/m or A/m. Without resumming other parts of the 
perturbation series, we cannot determine in which parameter range these terms can be neglected. For the following estimates we 
shall simply assume that we are in that range. 

The interpretation of gr as the physical interaction parameter is reinforced by the fact that the expression we have found ( I58> . 
depends on the detuning parameter A, the effective plasma frequency 77, and the atomic number density p in exactly the same 
way as the non-linear Kerr susceptibility coefficient Xg'^ for a dilute gas of two-level atoms, as obtained by Grischkowsky IitIi . 

with /ie as the dipole matrix element (for one component of the dipole vector) connecting the two states of the two-level atom. 
Rewritten in our notation, 

2 1^ '^^ 9r 
/^e = in— = -Ti ^ (60) 



which gives 



To compare the expressions, we write the 3D susceptibility, extracted from our effective action il6\ . in terms of the dimensionless 
bare coupling constant g, 

^(3)^1^^(^X. (62) 
16 pm^A'^ 3 
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If we assume that the 3D susceptibility renormalizes in the same way as the 2D dimensionless g (the renormalization comes 
from high k where the dimensional reduction from 3D to 2D is not relevant), then 

X?^^^^ = f (63) 

and since m « luq, our expression for x^'^-' is very close to Grischkowsky's. 

That the factor S/Stt « .85 is very close to one is of no significance, since our coefficient depends on the details of the 
ultraviolet cutoff. What is relevant, however, is that the two quite different approaches give essentially the same result for the 
susceptibility, and also that our result is independent of the bare non-linear coupling. 



6 SCALES, BOSE-EINSTEIN CONDENSATION AND TWO PHOTON BOUND STATES 

Based on the discussion of the effective photon-photon interaction in the previous sections, we will now consider some of 
the physical aspects of the formation of a two-dimensional photon fluid. We first summarize the important parameters that 
characterize the photon system, and then discuss the conditions under which two particularly interesting phenomena could 
occur: the formation of a two dimensional Bose-Einstein condensate, and the formation of two-photon bound states. We should 
stress that both these effects are essentially quantum mechanical, and cannot be described by classical non-linear optics. 



Important scales 

The strength of the mixing between the photons and the oscillators is given by the effective plasma frequency rj defined by 

(64) 

and mixing becomes important when rj^ ~ A^. From now on we shall restore factors of c and 7i in the formulas. 
The (unrenormalized) interaction strength is given by the dimensionless coupling constant g, 



1 i UK (Tj 
^ " (27r)2 16^ 



(1) ■ 



The importance of the non-linear loop corrections to the interaction strength is, for momenta ^ 2hinA, determined by the 
dimensionless parameter a, 



37r ;i^A 

(X = 

64 c^m^A 

and for large a the effective coupling constant is given by 

1 cm 7] 

Finally, for very small momenta, with 



(66) 



(67) 



fc2 

.gin— — ~l + a, (68) 
2mA 



the logarithmic term will become important and contribute to the renormalization of the interaction. 



Photonic Bose-Einstein condensates 



As mentioned in the introduction, one of the most exciting aspects of forming a photon gas, would be the possibility to study 
phase transitions. First, recall that there is no Bose-Einstein condensation (BEC) in a free two dimensional gas. However, the 
situation is different when the gas is in a trapping potential, where the transition temperature is given by fc^Tc ~ ^/NHQ f 18|, 
where N is the total number of bosons and ft is the frequency of the (harmonic) trapping potential. In an interacting gas the 
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situation is more complicated. The Mermin-Wagner theorem fl^ rules out true long-range order, but there is still the possibility 
of a Kosterlitz-Thouless (KT) transition, with the formation of a "quasicondensate" at a critical temperature ksTc ~ p2Df^^ /"m-, 
with the 2D boson density and m the boson mass. In the context of two-dimensional atomic BEC, the phase transition in 
a quasi two dimensional Bose gas with repulsive delta function interactions was recently analyzed by Petrov, Holzmann and 
Shlyapnikov | l^J. They find that in a trapping potential a true condensate will form well below Tc, whereas for intermediate 
temperatures T < Tc it will change to a quasicondensate with a spatially fluctuating phase. 

For a two-dimensional (dilute) photon gas a similar analysis should be relevant. We will make some comments on this in 
the discussion of a photon gas interacting with Rydberg atoms in the next section. For this discussion the expressions for 
the effective interaction parameter, i67\ and the two-dimensional cross-section i29\ are important. There are also important 
questions concerning the thermalization time and the possibility of regulating the effective temperature of the photon gas. With 
the parameter values discussed below the critical temperature is typically very much higher than the excitation energies. This 
indicates that a thermalized photon fluid would tend to form a Bose condensate rather than a normal fluid, and that the phase 
transition temperature cannot easily be reached by varying the parameter values of the photon fluid. 

A two-dimensional two photon bound state 

The scattering amplitude i32i always ha s a p ole for complex k on the second Riemann sheet, at = —2mEi,. This corre- 
sponds to a bound state, a "diphoton state" 12011 . with binding energy Eb given by 

2 

Eb = = |A|e^ . (69) 

Zm 

Somewhat surprisingly it looks like one could have a bound state for either sign of g^, and this is indeed the case for a fundamental 
delta function interaction |12|. We note, however, that for repulsive interactions (gr > 0), the bound state would occur at large 
k, far outside the range of validity of our effective theory. For attractive interaction (gr < 0) however, the presence of such a 
two photon bound state is a bona fide quantum mechanical effect with no obvious counterpart as a solution of the non-linear 
Schrodinger equation. 

When \a\ < 1, attractive interaction corresponds to negative A, and diphoton states could in principle form for red detuning, 
where the photons are stable against decay to the lower energy dipole mode. However, when a >> 1, which seems most relevant 
for physical realizations (se discussion below), weakly bound photon pairs would occur only for blue detuning. In this case we 
do not expect the formation of stable bound pairs due to the presence of the decay channel to the dipole mode. 

7 THE 2D PHOTON FLUID: PHYSICAL REALIZATIONS 

In this section we shall discuss two possible experimental scenarios for the formation of a 2D photon fluid. First, we consider 
microwave photons in a cavity filled with Rydberg atoms. This scenario is quite close to the one described in the previous 
sections, since the highly excited atoms are well described as harmonic oscillators with a small anharmonicity. As a second 
example we shall consider optical photons in a cavity filled with alkali atoms in their ground states. 

The question we shall address is the following: Can the parameters of these systems be tuned in such a way that there 
results a sufficiently large effective photon-photon interaction — mediated by virtual transitions within the atoms — so that a two- 
dimensional photon fluid can form within a Fabry-Perot resonator? One might expect that such a fluid arises after many effective 
photon-photon collisions within the cavity. Based on the analysis of the previous sections we reach an answer to this question 
which is qualified affirmative. The conditions for a photon fluid in thermal equilibrium to form may be met, but only under 
conditions where the parameters can be tuned to optimal values. For the Rydberg atoms these conditions seem more difficult to 
satisfy than for the optical photons. 

However, the conclusions concerning the formation of a photon fluid are based on simple order of magnitude estimates of the 
physical parameters. One should note that the expressions used in these estimates were found above assuming certain conditions 
(separation of physical scales) which are not necessarily met in the real system. Also the expression found for the interaction 
parameter 5,. is based on the resummation of certain large contributions, while other higher order contributions are left out. For 
this reason we can present only a preliminary evaluation of the conditions for a 2D photon fluid to form. A detailed evaluation 
of these conditions is outside the scope of the present paper 
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Rydberg atoms in a microwave Fabry-Perot cavity: Numerical estimates 



We first consider the case where microwave photons interact via Rydberg atoms excited near resonance, with both photons 
and atoms injected into a microwave version of a Fabry-Perot cavity. Consider a ribbon-shaped beam of Rydberg atoms passing 
through the gap between a pair of parallel, conducting plates, separated by a spacing Ll ~ Ao/2. The plates are partially 
transparent to microwaves, for example, by having a regular square array of small holes, and can thus act as the (highly reflective) 
mirrors of a microwave Fabry-Perot cavity. Microwave photons can be injected into the cavity from the outside, and their 
transverse distribution, after many photon-photon interactions mediated by the atoms, can be monitored in transmission. The 
resonance frequency of the cavity can be tuned to near a resonance frequency of the Rydberg atoms. 

For simplicity, let the plates be square in shape, with a typical transverse dimension of Lt » Aq, where Aq is the microwave 
resonance wavelength of the Rydberg atom transition of interest. The longitudinal mode number of the cavity is restricted to that 
of the fundamental longitudinal mode of the cavity (by our choice of spacing Ll = L — Xq/2 between the plates), but there can 
exist many possible transverse modes of the cavity, which are closely spaced near the fundamental longitudinal mode. 

For concreteness, we shall use as a guide the parameters of the experiment of the Paris group |21], where Rubidium atoms 
were excited to n = 50, where n is the principal quantum number of the Rydberg atom. The Rydberg atoms can be put into the 
"circular" state \n,l = n — l,m = n — 1), which have a maximal electric dipole moment, and thus a maximal coupling to the 
microwave photons. Since the typical size of such a Rydberg atom is given by the radii of the circular Bohr orbits a„ — ri^ao, 
where ag w 0.5 A is the Bohr radius, the atom can be quite large in size, e.g., a„=5o ~ 0.125 microns. The electric dipole 
transition matrix element for the ?T. = 50-^n = 51 transition is also quite large, and therefore the photon-photon coupling 
mediated by such atoms virtually excited near resonance should be correspondingly large. 

In the regime of high n, the energy levels of Rydberg atoms are almost equally spaced, 

Ry 2Ry ^ 3Ry , , 

where Ry — 13.6eV is the Rydberg constant. The dominant photon-induced transitions will be between the circular states. 
In, n — 1, n — 1) ^ In ± 1, n ± 1 — 1, n ± 1 — 1), where the wave function of the state In, n — 1, n — 1) is 



{2k) 



Sttti 



A— I -e"'^'' (2Kr)"-i (sin6l)"~i e*("-i)'^ , 



with K ~ 1 /nag. This system is very well suited to be modeled by an anharmonic oscillator model of the type discussed earlier. 
We shall determine the anharmonic parameter A, the oscillator frequency uj, and the oscillator mass M of this model in terms of 
the Rydberg atom parameters, by matching the both the energy spectra to terms quadratic in n and the dipole matrix elements 
between the levels n and n + 1. 

One might object to this procedure since in our previous calculations we assumed the oscillators to be in their ground state 
rather than in a highly excited state. In particular this means that the effect of virtual transitions to lower Rydberg states that 
would be present in Rydberg atoms with the SDectrum (l70> are not taken into account. However, it should be kept in mind that 
the detuning for transitions of Rydberg atoms to lower states will be substantially larger than those to the higher states. In any 
case, for the simple numerical estimates made in this section we shall ignore this effect, and assume that only upward transitions 
n — > n + 1 are important. This approximation will at most give rise to a numerical factor of order of unity in the estimate of 
the model parameters. For the same reason we shall neglect the fact that the oscillators of our model are three-dimensional and 
simply match the dipole matrix elements of the Rydberg states with those of the corresponding one-dimensional anharmonic 
oscillator. 

With the above matching procedure for "circular" Rydberg atoms prepared in a state with the initial principal quantum number 
n, we get the following result for the anharmonicity parameter A of the oscillator model 



K = --^ = ii. (71) 

However, this expression gives for the magnitude of the renormalization constant, |a| ~ ujo/{nA), and a therefore becomes very 
large, since we shall operate the system with very small detunings near resonance, so that A << ljq. The detuning A, however, 
must be much larger than the natural line width of the atoms |22], and in this limit we can use the A independent expression 
(I67> for the renormalized coupling constant. Note that in this limit the sign of the physical interaction is determined solely by 
the detuning parameter A - repulsive for red detunings A > 0, attractive for blue detunings, A < 0. Matching the dipole matrix 
elements effectively amounts to introducing a large oscillator strength, fosc = n^(2n — — 1)'^ in the expression for the 

effective plasma frequency rj — (q-^p/M)^^'^ — ^ [hce^ foscP/rneY/"^ = c{aefoscXcPatomY^^, where ae is the fine structure 
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constant, me is the electron mass, and Ac = 2ttXc is the Compton wavelength of the electron. Although, strictly speaking, the 
detuning A should be much greater than the effective plasma frequency ij, for the purposes of numerical estimates, we shall set 
A equal to 77. Also, we shall approximate mc^ w where wq is the atomic resonance frequency of the transition in question. 
With this we get, 



where Aq is the atomic resonance wavelength of the transition in question. 

As a first estimate of the size of Qr, consider the case of the the n = 50^n+l = 51 "circular" Rydberg transition, where 
mc^/h « Wo = 27r X 51 GHz, with p ~ 1000 atoms per cm^ and fosc = 105. This yields an effective plasma frequency of 
7] ~ 5.15 X 10^s~^. For this transition, assuming that A ~ 77, we find that Qr 6.32 x 10"'^. 



Cavity dimensions 

To fit the frequency of the Rydberg transition n = 50 — > 51, i.e., loq = 51 GHz, we assume the dimensions of the cavity to 
be chosen as follows, 

Ll = ^Aq ~ 3 mm and Lx ~ 100 cm, (73) 

with Ll = L as the longitudinal extension of the cavity and Lt as the extension in the transverse directions. The limits for the 
transverse momentum kr = k are therefore 

^ <kT< (2toA)^/^ . (74) 

For a detuning A ~ 27r x 2 MHz, this implies that there are around 10 nonrelativistic transverse modes in the cavity which 
can interact via the Rydberg atoms. Higher modes will be available and they will also interact via the atoms, but the energy- 
momentum relation for these higher modes will be modified relative to the non-relativistic expression. However, even with the 
limited number of modes available one should in principle be able to see the formation of a photon fluid in a Bose condensed 
state. 



Tiie effective piioton-piioton collition rate and tiiermalization 



We showed earlier that the total elastic cross-section in 2D is given by 

,2 



a^' « f^. (75) 
The "reaction rate" for photon-photon collisions in 2D is therefore 

cr"'i^T = a^' = —gl. (76) 

m m 

The collision frequency WcoH for photon-photon collisions in 2D is 

WcoH = P2D(t''^Vt = P3dLl(T^''Vt , (77) 

where p2D is the 2D photon number density, and psn is the 3D photon number density. Now the number of photon-photon 
collisions N that occur within a cavity ring-down time r is given by = lUcoUT, where, given the quality factor Q of the cavity, 
T = Q/coq. Solving for the number density of photons pso needed for Ncoii photon-photon colUsions, one obtains 



P3D = 87r2-^Ao^ (78) 

Q9r 
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We note that due to the smallness of the interaction parameter (gr ~ 6.3 x 10 there is a large numerical factor in these 
expressions. If we assume ~ 10 and a rather high quality factor of Q — 10^, 

-— = 2.5 X 106 (79) 

Q9r 

and the corresponding value of the photon density therefore gives a KT transition temperature far above the low-momentum 
regime of the interacting photons, 

ksTc « w 2 X 10*^ hujo (80) 

m 

corresponding to a temperature Tc — 4 x 10^ K. Thus, the low-energy photons are typically in a temperature interval where a 
Bose condensate will form, rather than a normal fluid. 
The microwave intensity inside the cavity is 

Iin=P3DchujQ (81) 

and with the given parameter values we find I^o w 1.1 mW cm^^. To convert from the inside-cavity to the outside-cavity 
intensity, we use the relation lout = I-in/Q- Note that the incident power lout needed in order to form the photon fluid scales 
inversely as the square of the quality factor, i.e., lout ~ Q^^, since /i„ is inversely proportional to Q. 

Based on the above, we estimate that about 1 . 1 microwatts of incident microwave power is needed in order to get around ten 
photon-photon collisions within a cavity ring-down time of 30 microseconds. Under such circumstances, we expect that a 2D 
photon fluid should form inside the cavity. 

The value of Q ~ 10^ is feasible (Q ~ 10^ by means of superconducting cavities have already been achieved lEill '). Note, 
however, that if A > the situation rapidly deteriorates due to the dependence on the factor (rj/ A)^ of g^. 



Chemical potential and speed of sound of ttie photon fluid 

From the previous section it might seem that it is always possible to form a photon fluid just by increasing the intensity of the 
incident microwaves. This conclusion is however not justified, since a too high intensity will make the fluid so dense, and the 
mean kinetic energy is so large, that the non-relativistic approximations are no longer valid. To get a rough estimate of when 
this happens, we use an effective Ginzburg-Landau (GL) type theory that should be applicable when the fluid is very dense. The 
relevant potential is given by 

V^-ti\cj)\^ + ngr—\^f^O. (82) 
2m 

where we introduced a chemical potential which, in the Thomas-Fermi approximation, is related to the photon density by the 
relation 

t^^'^9r7^P2D, (83) 

2 m 

which is obtained by minimizing V with respect to |0p — p2D- The GL theory defined by V implies sound waves in the photon 
fluid with a speed given by 



1/2 

Vs = [-) ■ (84) 

\m/ 



We can also estimate the mean velocity of the photons using the virial theorem, 

{Ekin) — ^ P2D — {Epot) — PP2D (85) 

which gives V gamma = V^V^. 

Imposing the consistency condition that Vs — (/i/m)^/^ < c, we obtain an upper limit on the 2D photon number density that 
P2D < 2m^c^ /h^Trgr- Combining this with the formation condition that iV >> 1, and using lvq — inc^/h, we obtain the upper 
and lower bounds on p2D 

m^c^ 2m^c^ 

-J— « P2D < -2 . (86) 

Ti Qg^ n irgr 
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We therefore conclude that a self-consistency condition for the formation of the photon fluid is 

Q» ^gr' ■ (87) 

Thus, the choice of a high-Q resonator is an important criterion for the formation of a photon fluid. For the case of the Rydberg 
transition n — 50 ^ 51, we see that it would be necessary for Q >> 10^, so that the choice of Q ~ 10^ may be sufficient. 

Alkali atoms in an optical Fabry-Perot cavity 

The apparent advantages of using microwaves and Rydberg atoms, rather than a more standard setup with lasers and e.g. 
Rubidium atoms, is twofold. First the cavity dimensions are not too small, and second, the large dipole moment of the atoms 
implies a large nonlinearity. The obvious disadvantage is that the Rydberg atoms are much harder to create and manipulate. 
Also since we have shown that the anharmonic parameter A drops out of the renormalized coupling constant for small A, it 
is of interest to consider an optical Fabry-Perot cavity with alkali atoms in their ground states, for which there is only a single 
transition to one excited state with a very strong oscillator strength of the order of unity, which almost completely exhausts the 
/-oscillator sum rule. For alkali atoms, a two-level model for the atom is therefore a more appropriate one 1 17]. However, 
the renormalized coupling constant for photon-photon interactions, and the resulting rates for collisions found in the previous 
section, should still apply here, provided that we remember to use the condition that the detuning A be comparable to the 
effective plasma frequency rj for the expressions for collision rates. 

Consider the case of a fundamental longitudinal mode Fabry-Perot cavity with Ll ~ Aq /2. In this case, the Q of the cavity is 
also approximately its finesse. High quality mirrors with Q ~ 10^ are commercially available 1 24]. For a strongly allowed alkali 
transition, fasc — 1- If we assume a density Patom — 10^ rubidium atoms per cm'^ (rubidium atoms have a transition wavelength 
of Aq = 780 nm), this leads to an effective plasma frequency 77 ~ 27r x 200 MHz (which is comparable to the Doppler width of 
rubidium atoms at room temperature). With A ~ 77 this gives ~ 12 x 10^"^. The minimum required Q is therefore around 
130, which is easily satisfied. Thus, if we set TV ~ 10 and Q ~ 10^, from Ea.( l81> we obtain I^u ~ 10^ W cm^^, or an outside 
power requirement of a 0.1 Watt per square millimeter incident on the Fabry-Perot, which is feasible. 

One consequence of going from microwave to visible wavelengths, is that the number of nonrelativistic transverse modes 
given byEq.lllll can be much larger For the given parameter values one finds that there are around 200 nonrelativistic transverse 
modes which can be coupled together. Also the estimated KT transition temperature is relatively lower, with fcsT'c w 0.7 fiujo. 

High-Q microspheres of glass and the formation of a 2D photon fluid 

Let us finally comment on a third possible way in which a 2D photon fluid could form. Extremely high-Q optical cavities with 
Q ~ 10^" have been fabricated out of low-loss, small-diameter glass microspheres j^l- A 120-/im diameter glass microsphere 
immersed in superfluid helium has been observed to exhibit a nonlinear, dispersive bistable behavior with a threshold power 
of as low as 10 microwatts, due to the intrinsic Kerr nonlinearity of the glass which constitutes the microsphere 1.26.] . Due 
to the curvature of the microspheres the photons are tightly confined to propagate within an optical wavelength or so of the 
two-dimensional spherical surface, i.e., within its "whispering-gallery" modes. This tight, two-dimensional confinement of the 
light should also allow an effective dimensional reduction from 3D to 2D in the photon degrees of freedom, which differs from 
the mechanism of dimensional reduction in a planar Fabry-Perot cavity discussed here in this paper Thus a two-dimensional 
photon fluid could also in principle form just within the inside surface of the microsphere, due to the photon-photon scatterings 
mediated by the intrinsic Kerr nonlinearity of the glass, and may have already been indirectly observed in the experiments of 
It should be noted that the change of topology from that of a plane to that of a sphere should not make the formation of the 
2D photon fluid impossible in principle. 

8 CONCLUDING REMARKS 

In this paper we used a microscopic approach to study the effective interaction between photons confined to a quasi-two- 
dimensional cavity filled with a non-linear medium. With the atoms modelled by a collection of non-linear Lorentz oscillators, 
the interaction was studied to second order in perturbation theory. The linear problem was first solved by decoupling the 
"polariton modes", and we described efficient ways to do this, both by use of a path integral approach and by use of exact 
diagonalization, as shown in Appendix A. A main motivation was to examine the description of the effective interaction, induced 
by the non-linearity, as a 2D short range (delta function) interaction and to determine the renormahzation scale associated with 
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this interaction. We noticed several interesting complications. For red detuning relative to the oscillator frequency, the expression 
found to second order fits with the form of a normalized delta-function potential. However, the renormalization parameter is 
typically large and a resummation of "dangerous" terms has been performed in order to extract an effective, renormalized 
interaction strength. This way of including higher order contributions raises questions of whether other contributions, not 
included here (since they are parametric ally small), may also be of importance. However, additional support for the expression 
found for the renormalized interaction parameter comes from other evaluations of the non-linear susceptibility of a gas of alkali 
atoms. For blue detuning, we discussed the possible formation of photonic bound states, and noticed the possibility that photons 
might "decay" into oscillator excitations. 

The last part of the paper was devoted to a discussion of possible experimental scenarios, where a photon fluid may form 
in a Fabry-Perot resonator In these cases, two-dimensionality is obtained by constraining the photons in one direction to their 
fundamental mode. By means of order of magnitude estimates we examined the constraints on the physical parameters in order 
for a photon fluid to form. There is a potential conflict between the need to have a high density of photons in order to obtain a 
satisfactory collision frequency and the possibility of having a too high density measured in units of the photon Compton wave 
length, since the effective photon mass is typically very small. For photons in a microwave cavity interacting with Rydberg 
atoms the corresponding constraints are serious, although by choosing optimal parameter values the conditions for creating a 2D 
photon fluid can probably be met. Such an experimental realization is otherwise attractive, as the model of photons interacting 
with Rydberg atoms is closely related to the oscillator model used in the paper. For photons in an optical cavity, the constraints 
on the physical parameters seem less severe, although the use of such small scale cavities may pose a larger experimental 
challenge. We noticed that in both cases the typical temperature associated with a Kosterlitz-Thouless phase transition lies well 
above the energy scale associated with non-relativistic 2D photons, and for such low energies a condensate will typically form 
if the collision frequency is sufficiently large on the scale set by the ringdown time of the cavity. In addition to the two possible 
experimental realizations discussed in some detail, we also suggest that other types of realizations of the 2D photon fluid may be 
possible, and in particular we pointed at the use of high-Q microspheres of glass as being potentially interesting in this context. 

Let us finally stress the point that in this study we have made some simplifying assumptions, in particular about a clear 
separation of the physical (energy) scales associated with the formation of a 2D photon fluid. In real systems, these conditions 
may not be well satisfied, and in the discussion of microwave photons interacting with Rydberg atoms we noticed that one may 
have to tune the frequency so close to the resonance value that the photons become strongly mixed with the oscillator degrees 
of freedom, i.e., they are genuine polaritons rather than photons. For this reason we consider the order of magnitude estimates 
applied to these realizations only as preliminary ones. A more detailed study is needed to settle more firmly the conditions for 
the creation of 2D photon fluids in such small Fabry-Perot cavities. 
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APPENDIX A . ALTERNATIVE DERIVATION OF THE EFFECTIVE PHOTON ACTION 

In this Appendix we derive the effective action il4\ and (I16> by an alternative method which does not require any field 
transformation of the type (I13> which is non local in time. In addition to verifying the soundness of the simple procedure used 
in the text, the present derivation explicitly shows that the effective action is the sum of the two contributions in J14> and (I16> 
corresponding to the two frequencies uj±. 

The most natural object to consider given the action Q, would be the generating functional Z[J^, J^], which can be written 
as a Feynman path integral over the fields A and R. By a change of variables one can diagonalize the quadratic part of the 
action, and then use perturbation theory to calculate Greens functions. However, becuause of the time derivative coupling @, 
this will give rise to expressions which are nonlocal in time, since this procedure fails to correctly describe the normal modes. 
An obvious way to proceed is to switch to a Hamiltonian formalism, do an explicit diagonalization of the quadratic piece of the 
Hamiltonian by a symplectic transformation on the phase space, and finally express the nonlinear term in the new variables. This 
method is however cumbersome, and we now present a simpler approach based on the path integral. 

The essential step is to use an alternative form of the action which is related to Q by a Legendre transformation. Rather 
than just writing down this transformation, we shall go back to the derivation of the Feynman path integral, since this will 
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illuminate the physical meaning of the new variables. The usual Hamiltonian description of ([T) would comprise two canonical 

pairs (A, —E) and (_R, P), where P — MR + qA. Since we are dealing with a coupled system of matter and radiation, it is 
rather natural to introduce the electric displacement vector by 

D = E + qpR = \7xC, (88) 

where the last identity defines the "displacement vector potential" C provided that there are no macroscopic charges, i.e., 
V • Z? = 0. The variables D and C were used earlier by Hillery and Mlodinow in the context of quantization of non-linear 
electrodynamics iEtIi . It is easy to show that C forms a canonical pair with the magnetic field strength B. If we furthermore 
introduce the rescaled oscillator fields 

Q = ^/pMR (89) 

n = -^=R, 

VpM 

the canonical commutation relations take the form, 

[B,(r),Q(r')] = inS,,5^{r~^) (90) 
[Q,(r),n,(r')] = m^S^r-r'), 

and the corresponding quadratic Hamiltonian density is, 

^ i(^B^ +D^) + ^[li^ + iu^ + rj^)Q^] ~VD-Q, (91) 

where again D = V x C, and 77^ = q^p/M. We now consider B and 11 as the momenta and derive the path integral expression 
for the generating functional 



Z[JdJr] = j i?(7I?Qe"5['^''5]+«/(./i, i5+/«.i?) 



(92) 



by first writing a phase space path integral and then integrating over the momenta B and n. As we shall see below, the electric 
displacement field C will in the limit of weak coupling correspond to the photon, and the oscillator Q to the oscillating dipoles. 
At finite coupling these two modes mix and form two effective propagating "polaritron" modes. In Fourier space, the quadratic 
part of the Lagrangian takes the form 

where P^^ = 5ij — kikj is a transverse projector Since in this formulation, the dipole coupling term does not involve any 
time derivative, the quadratic part of the action can be diagonalized by a unitary transformaton which is local in time (although 
nonlocal in space): 



Ci\ ( cos,QkPf. -isinOkei.kkk 



QiJ \ismdkeijkkk cosOkP^j 




(94) 



where 



cos9k^\l^ ^ (95) 

with k'^ — P and, uj±{k'^) given by ( fT2t . 

The resulting action is most conveniently expressed in spherical components defined by R„i = ■ R etc. , with e* 
Smm', k ■ im — and m = ±1. We get 



2 

m— ±1 
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For the interaction term, we must express R in terms of A'-^^ using ( I94K 

R,„{k) = (cosekA<^\k)+smekAi^\k)) 



where we also redefined A^^^ (k) ±^ii^ (k). 
The final result for the generating function is. 



with S[A^+'> , v4(-)] = 5(2) + ^^''^ exphcitly given by 



£(2)(1(+),1(-)) 



1 f cPk f duj 



S. 



(4) 
eff 



2 y (27r)3 J 27r 



pXM^ 



n 

n=l 



(2^)4 



(27rf,5(^c.„),5(^4) 



X ^(fci, wi) • i?(fc2, W2) -R(/c3, W3) • -R(A:4, W4) , 



(97) 



(98) 

(99) 
(100) 



where i? is given by j97r Note that the interaction is non local becauce of the fc^ dependence in the denominators in the cos 0k 
and sin Ok factors. Using the identities 



cos Uk = 



sin^ 6k 



(101) 
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and substituting the expression for R we recover ([T^ . 



APPENDIX B. SCATTERING AMPLITUDE TO SECOND ORDER 



We perform here a more detailed evaluation of the second order contributions to the scattering amplitude. The contributions 
from low and high momenta of the intermediate states are examined, and we consider in both cases only the leading contributions 
in terms of the small quantities k/m. A/m, and 77/ A. The contributions come from diagram a) in Fig. 2, with the intermediate 
lines corresponding to either the photon or the dipole mode, with transverse momenta running from zero to the cutoff \/losc, 
where lose is the physical distance between the oscillators. There is also a sum over all longitudinal momenta, in terms of the 
discrete momentum variable n, and a sum over the polarization vectors of the fields. 

The full expression for the interaction Hamiltonian is 



hm = n 
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+ "3,712 +"4 



ni+n2+n3,n4 



X ^ n (^'■' t) • (fe, n2, t)^fe) (^3, «3, t) ■ E^'"^ (ki, n4, t) 



where K = {2Tr)^ Arj^/ {8 pL) and A± (fc, n) = [Lu±{k, 



0) 



(102) 



Tlrir 



1, .., 4) is the summation variable for 



the longitudinal mode and ir for the photon/dipole mode. From this expression the T-matrix and the corresponding scattering 
amplitude can be determined. The T-matrix is related to the scattering matrix by S*/,; = Sfi~ 2TTi6{Ef — Ei)Tfi, where Ef and 
Ei are the energies of the final and initial states, and to second order Tfi has the form 



Tf, = {k[k'^\H,nt\kik2) + \j ^'Pi^'P2^ 



{k[k'2\Hint \PlP2){PlP2\Hint\klk2) 



) + uj{k2) - Lo{pi) - uj{p2) + i e 



(103) 



where we here have simplified the notation by suppressing the sum over field modes, longitudinal momenta and polarization 
variables in the intermediate state. Galilean invariance implies the following form 



Tj,^5{Kj~K{)T{ks,h) 



(104) 
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where Kf and Ki are the sum of momenta for the outgoing and incoming particles and kf and ki are the relative momenta. For 
pure s-wave scattering (delta function interaction) T only depends on the magnitude of the relative momentum, 

T{kf,h)^T{k^) (105) 

and is related to the s-wave scattering amplitude through 



/ = -(2^)^my^r(fc2) (106) 
For the present case, the first order contribution to T is 

^ 27rm 16(27r)3mpA4 

With the interaction matrix elements evaluated and the sum over polarization performed, the expression for the second order 
contribution to T is 



niri2 1112 



X / d^p 



{{Oji, {p, ni))2 - a;2)2 + 772^^2 ((^.^ n2))2 - ^2)2 + tj2^2 



X (1-1 . P\ .Vl ' \ 



2p^+n\ml) \ 2p'^+nlmlJ 2tj_(fc, 1) - a;ii(p, ni) - ^^^^(p, 712) + 



(108) 



with C ~ K/ (4(27r)''mA2). Here the reference frame is chosen where the total transverse momentum vanishes. The energy of 
both incoming and scattered photons is lo- (fc, 1), with the photons in the lowest longitudinal mode (and where red detuning has 
been assumed). Since these are low momentum photons, we have a;„ (fc, 1) ~ m + k"^ /m. 

A simplification may be introduced by considering the effect of mixing between photons and dipoles, described through the 
parameter 77. This mixing affects the expression in two ways. The first one is through the rj dependent term that is explicit in 
(I108> . This term is important when the particle energies are close to the resonance value wq- The other one is the indirect one 
which enters through the particle energy uj±. When summing over both the photon and dipole modes in the intermediate state 
this rj dependence is less important. It only affects the lowest longitudinal mode and only for transverse momenta in a small 
interval p « A. We will neglect this effect, which can be viewed as higher order in 7/2 and use the energy expressions for the 
uncoupled modes. 



^- = XlP^ 



= ujo (109) 

(Note, with this redefinition is the photon mode also when higher in energy than the dipole mode.) We consider now 
separately contributions from the different modes in the intermediate states. 



Two-dipole intermediate states 

Since w+ is independent of p there is no suppression of the integrand for high momenta and the integral over p and sum over n 
are divergent without the momentum cutoff. This cutoff we set to Pmax for the transverse momentum and N = PmaxL/ir for the 
longitudinal momentum. The connection with the oscillator density is p = l/l^^c — {Pmax / t^Y ■ Since we assume L » lose 
most of the contribution comes from large transverse and longitudinal momenta where the dimensional reduction no longer is 
effective . We note that the only momentum dependence now sits in the polarization factor, which we approximate by 

1 p2 1,3 

1 - - , , 9 =1-7: sin e^- (110) 

2 p2 + n^ml 2 4 

where 6 is the angle between the momentum vector and the 2;-axis. In the last term we have approximated sin^ 6* by 1/2 which 
is correct when we (for large momenta) can treat the longitudinal momentum as a continuous variable. 
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With these approximations the momentum sum and integral can trivially be performed, and the contribution is 

T2„ = -^C'nNp' ^ 



_ ■■'■r-max A 4- 2 
32 Ar/^CJQ 

We note that the sign of this ^-independent contribution depends on the sign of A, the detuning parameter. 

The photon-dipole contribution 

We now consider the terms where one of the intermediate modes is a photon mode and the other a dipole mode. The mo- 
mentum integral is convergent, with the main contribution coming from low momenta. For p « nm we approximate the 
polarization factors by 1, and T is simplified to 



dp 



(112) 



{p^ — uJq)"^ + tj'^lUq p + a — m. — ie 

The leading contribution (in A/m) comes from the term n = 1. Evaluation of the momentum integral gives for red detuning 
A > 0, 

9 1 A^jyS 



To 
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128 (27r)4 mp2A5 
8 27r pA ' 



-l^!^r, ("3) 



and for blue detuning 



_ 9 1 A?Tf 9 1 K^rf 



128 (27r)5 m/92A6 512 (27r)4 mp^A^ 

One should note the difference between red detuning (A > 0) and blue detuning (A < 0). In the latter case the integration path 
passes a pole which gives rise to an imaginary part. That is not the case for red detuning. 



Two-plioton intermediate states 

With two photons in the intermediate state the main contribution comes from the term m = 112 = 1 which has a pole at 
= k"^. The low momentum approximation, uj-{p'^) w m -|-p^/2to can be used, and the expression for T simplifies to 

9 f°° 1 1 

This gives as leading terms for red detuning 

^ 9 1 A^r^s ,^ . , 9 1 A^T/S 

^2c = 77^7:^^:7^ (log + 



1024 (27r)5 mp2A8 ' ^ 2mA ' 512 (27r)4 mp'^A^ 

3 1 Ary^ , k^ .,31 Ajj\ ^ 

' (log— — -Z7r) + — — Ti (116) 



64 (27r)2pA4' '=2mA ' 32 (27r) pA 



For blue detuning the expression is 



^ 9 1 A2,7« „ e ,91 h?ri 

^2c = /o_N.^^-9As (log;r:^-»7r)- 



1024 (27r)5 mp2A8 ' ^ 2m|A| ' 256 (27r)5 m/92A6 
3 1 Ar?" fc2 ,31 Ar;2~ 



^64 (27r)2 pA4' ''2m|A| ' 16 (27r) pA2 
The pole atp2 = k'^ gives as expected a term that depends logarithmically on fc2. 



(log - *'^) - 777 7^3^2 ^1 (117) 
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The full expression 



When adding the contributions we find for A > 0, 



_ , 3 1 A774 fc2 Stt A 9 1 A77 \ 

^ ~ '^64 (27r)2 pA4 ^ °S 2^ ~ '''' ~ 64 ~ 32 (2^^ ) ^ 

3 1 Ary^ , fc2 , 37r A \ ^ 

,64 (2^^^^°^ 2;;^ - - 64^^ ) ^''^^ 

where in the last expression we have assumed r/ << p/rn^ and have left out the small term. For A < we find 

^ / 3 1 Aif fc2 , 37r A 3 1, ,W\^ 

- (,64(2^^(^°S^-*^)-^;^ + Y6(2^(^+''')^j^^ 

3 1 A7y4 p _ 37r A . 3 1 A772 \ ^^^^^ 



64 (27r)2pA4' '=2to|A| ' 64 m2A 32 27r pA2 

We note that whereas the expression for red detuning agrees with the expected form for a renormalized delta function interaction, 
the expression for blue detuning contains an additional imaginary part. This is small compared with the constant real part, but 
has nevertheless some significance. As noted above it arises from the pole in the energy factor in the case where one photon and 
one dipole is present in the intermediate state. It is therefore related to the possibility of real scattering of the two photons into a 
photon and a dipole excitation, as discussed in the text. 
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